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P/7s in Metals Caused by Collision With Liquid Drops 

and Soft Metal Spheres 

Olive G. Engel 

An equation is developed to give pit depth as a function of collision velocity for pits 
formed in soft to medium-hard metal plates as a result of collision with liquid drops. The 
rear face of the target plate must be a free surface. The plate thickness must not be less 
than 1.5 to 2.0 times the drop diameter nor greater than 4 to 5 times the drop diameter. It 
is shown that, under the same conditions on the target plate, the equation is also applicable 
to pits formed in soft to medium hard metal plates as a result of collision with spheres of 
the same metal that flow radially as a result of the collision. Pit-depth- versus- velocity 
data obtained in other laboratories were used to test the equation. Metals used as targets 
were copper, 1100-0 aluminum, 2024-0 aluminum, lead, steel, soft iron, and zinc. Mercury 
was used as the drop liquid against copper, aluminum, lead, and steel. Water was used 
as the drop liquid against copper, aluminum, and lead. Spheres of copper, aluminum, 
lead, soft iron, and zinc were used against targets of the same materials, respectively. The 
equation can be used to calculate the dynamic compressive yield strength of soft to medium- 
hard metals. 



1. Introduction 

Current research on pits produced in high-speed 
collisions includes: (a) the investigation of solid 
projectiles impinging against solid targets (artillery 
experiments); (b) the investigation of solid projectiles 
impinging against liquid surfaces (water entry 
problems); (c) the investigation of liquid drops 
impinging against solid surfaces (high-speed rain- 
erosion research); and (d) the investigation of liquid 
drops colliding with liquids. The work that is 
described in this paper was initiated as part of a 
high-speed rain-erosion research program and the 
entire program was conducted under the sponsorship 
of the Materials Laboratory, Directorate of Labora- 
tories, Wright Air Development Center. 

It is difficult to test for the rain-erosion resistance 
of structural materials at very high collision velocities 
because of the problems involved in accelerating 
either waterdrops or test specimens to the velocities 
in question. It has been suggested that it may be 
possible to bypass these problems by using drops of 
high-density liquids instead of drops of water for 
the erosion tests. 1 To develop this idea into a 
reliable test procedure, it is necessary to know the 
corresponding velocities-for-equal-damage when the 
test specimen collides with a drop of a high-density 
liquid and when it collides with a drop of water. 
The pit-depth-versus-velocity equation presented 
in this paper was developed to provide this 
information. 

It has also been found that this pit-depth-versus- 
velocity equation is applicable without change to 
collisions of spheres of the soft, ductile metals with 
targets of the same metal in those cases in which 
the sphere flows as a result of collision with the metal 
plate when the collision velocity is as high as 5,000 
ft/sec. Pellets of the soft, ductile metals appear, 

1 This idea was suggested by the rain-erosion research group working at Convair 
Division of General Dynamics Corp. in San Diego, Calif. 



therefore, to behave as though they were liquid drops 
when they impinge against a solid surface at this 
velocity. 

The experimental work reported in this paper was 
done in other laboratories. 2 :i 

2. Collisions With Liquid Drops 

The damage done to solid materials that collide 
at high speed with liquid drops depends both on 
the properties of a drop of liquid in collision and on 
the characteristic properties of the solid. 

2.1. Damaging Properties of a Liquid Drop in High- 
Speed Collisions 

As the relative collision velocity changes, the 
effects of collisions with a liquid drop also change. 
In high-speed collisions with the planar surface of a 
solid the liquid drop acts as though it were a hard 
solid sphere, but, unlike a. sphere of hard solid 
material, it undergoes an ultrarapid radial flow 
outward about the point of impingement [1,2]. 4 
These damaging properties of an impinging liquid 
drop vary in intensity depending on the density of 
the liquid of the drop, on the relative collision 
velocity, and on the extent to which the solid surface 
yields under the blow. At any given impingement 
velocity, collision of a ductile metal plate with a 
liquid drop does not produce as deep a pit as collision 
of the same metal plate with a hard solid sphere 
would produce. This is because when a planar 
solid collides with a liquid drop part of the collision 
energy is transformed into the radial flow of the 
liquid. 

2 The pit-depth-versus-velocity data reported for collisions of solid targets 
with liquid drops were obtained by the rain-erosion research group working at 
Convair Division of General Dynamics Corp. 

3 The pit-depth-versus-velocity data for collisions of spheres of the soft metals 
that flow as a result of the collision were obtained by Partridge, VanFleet, and 
Whited at the University of Utah (fee reference [11].) 

4 Figures in brackets indicate the literature references at the end of this paper. 



229 



2.2. Response of Solid Target Materials 

The extent and type of damage produced on solids 
as a result of collision with liquid drops depends 
strongly on the characteristic properties of the 
solids [1]. If the damage is to be described mathe- 
matically, it will be necessary to develop separate 
equations for the damage produced on solid mate- 
rials that have widely different properties. 

2.3. Dimensional Analysis of the Damage 

To develop an equation that will give the depth 
of pits produced in high-speed collision between the 
planar surface of a solid and a liquid drop is a 
difficult problem from a theoretical standpoint. 
Dimensional analysis has proved to be a useful 
tool in the solution of difficult problems in the past. 
It is used in this paper to develop a pit-depth- 
versus-velocity equation. 

The method of dimensional analysis has been 
discussed by Buckingham [3, 4], Bridgman [5], 
Birkhoff [6], and others. If a physical process can 
be described by physical quantities of n different 
kinds, and if none of the quantities involved in the 
process has been overlooked, then the process is 
described by the equation: 



/(Qi, &, 



Qn, r l, T 2<> 



)=o, 



(1) 



where the Q's are the physical quantities involved 
in the process, and the r's are ratios. The Q's are 
quantities of different kinds. If several quantities 
of any one kind are involved in the process, they are 
specified by the value of any one of them, and by 
ratios of the others to this one. These ratios are 
the r's of eq (1). 

A certain number k of the Q's are selected as 
fundamental. The k selected Q's comprise a possible 
set of fundamental dimensions and the remaining 
Q-quantities can be expressed in terms of them. 
After the k fundamental Q's, now designated as 
S's, have been selected, the remaining Q's are 
designated as P-quantities. In a mechanical system 
the number of the ^-quantities is three because the 
total number of dimensions required to express any 
of the Q-quantities involved in a mechanical process 
is three. These dimensions are usually mass, length, 
and time. 

Equation (1) may be put in the form 
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where the it's are independent products of the argu- 
ments Q, and are dimensionless in the fundamental 
units, or in the form 

ir 1 =F / (tt 2 , 7t 3j . . . 7r n - k , r u r 2 , . . . ) (3) 

in which wi is any one of the t's. 

A 7r-product is formed by multiplying each P- 
quantity by the S-quantities, which are raised to 



whatever powers a, 0, ... it are required so that 
their units will just cancel those of the P-quantit}^ 
to make the 7r-product dimensionless. That is, 



*r=(Q?-$ 



Ql)P. 



(4) 



From a consideration of the equation for the 
pressure that results when a solid surface runs into 
a liquid drop at high speed [2], of the expression for 
the radius of flow of the liquid about the central 
point of the collision [2], and of the well known 
equation for the shear stress between layers of 
liquid in laminar flow, the quantities (see table 1) 
that should prove to be important in determining 
the damage that results from the collision of a solid 
surface with a liquid drop are: c, c\, p, p ', V, pi, y, y' ', 
d, and 8' where c is the speed of sound, p is the dens- 
ity, V is the relative impingement velocity, p is the 
viscosity, 71s the surface tension, d is the diameter of 
the liquid drop, and 5' is the damage parameter. 
Primed quantities refer to the solid material; 
unprimed quantities refer to the liquid of the drop. 
In addition to these quantities there is the energy 
per unit volume put in, E' , and the energy per unit 
volume returned, e' , during the collision. E f is the 
energy per unit volume that the solid material can 
absorb without nonrecoverable deformation or frac- 
ture; e r is the energy per unit volume that the solid 
material can return. The ratio of these quantities, 
e' 'jE' , is a measure of the resilience of the solid 
material. 

Table 1. The three fundamental S-quantities, the P-quantities, 
and the r-ratios involved in rain erosion damage 



Quantity 


Dimen- 
sions 


S 


P 


r 


Dimensionless 
products, 7T 


Impact velocity, F..._ 
Drop diameter, d.. .__ 
Density of liquid, p - . _ 

Viscosity of liquid, ju. - 

Surface tension of 
liquid, 7 

Velocity of sound in 
liquid, c 

Damage parameter, 5'. 

Yield or rupture en- 
ergy density of the 
solid, E' 

Ratio of sound veloc- 
ities, c/c' 

Ratio of densities, 
pIp' 

Ratio of surface ten- 
sion and interfacial 
tension, 7/7' 

Resilience, e'/E' ._ 


LIT 
L 

M/L* 

MILT 

MIT* 

LIT 

L 
MILT"- 


Si 


Pi 

Pi 

Pz 

P 4 
Ps 


n 


7Tl- 1= prfV/jU 

Reynolds 
Number 

TT2- 1= pdV 2 ly 

(Weber 

Number") 2 
tt3=c/I 7 (Mach 

Number, 

liquid)- 1 
Ti=8'/d 

T & = E'f(pV2) 





























Of these quantities V, d, and p are arbitrarily 
selected to be the ^-quantities, S ly S 2 , and S 3 , 
respectively. Considering both these S-quantities 
and the remaining quantities, there are four pairs 
of quantities of the same kind: c, c' '; p, p \ y, y r ; 
and e', E'. These are expressed as four ratios, 
r u r 2, 7*3, and r 4 , respectively. One member of each 
of these pairs is a P-quantity if it has not already 
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been chosen to be an ^-quantity. Also any other 
of the tabulated quantities that has neither been 
expressed in a ratio nor chosen to be an S-quantity 
is a P-quantity. The P-quantities are then: 

fx, y, c, 8', and E' (see table 1). 

The 7r-products are formed by multiplying the 
P-quantities by the product of powers of the S- 
quantities. Each S-quantity is raised to whatever 
power is required to make the x-product dimension- 
less. For example, the viscosity fi has the dimen- 
sions M/(LT) where M is mass, L is length, and T 
is time, and, therefore: 

tc{- 1 =p d Vlfx=(M/U) (L) (L/T) (LT/M). (5) 

Surface tension y has the dimensions M/T 2 and, 
therefore, 

t 2 -' = p d V 2 /y=(M/L*) (L) (L 2 jT 2 ) (T 2 /M). (6) 

7r 3 , 7r 4 , and 7T5 are formed by a similar procedure. 
7Ti _1 and 7r 2 _1 are the dimensionless Reynolds Number 
and the square of the Weber Number, respectively. 
7T3 -1 is a kind of Mach Number that gives the ratio 
of the relative impact velocity between the liquid 
drop and the solid surface to the speed of sound in 
the liquid of which the drop is composed. 

The damage parameter <5' may be chosen to be a 
length, a volume, or an area. A length is the easiest 
quantity to measure experimentally. Taking <5' to 
be a length, eq (3) for the damage caused by high- 
speed collision of a solid surface with a liquid drop is 

l'ld=F'[pdV/v, P dV*/y, V/c, E'UpV 2 ), c'/c, 

p'/p, 7 77, e'lE']. (7) 

It was pointed out in section 2.2 that it will be 
necessary to develop separate equations to describe 
the damage produced in liquid-drop-solid-surface 
collisions for materials that have widely different 
properties. The damage produced on nonrubbery 
materials that fail by plastic yield but without frac- 
ture is the easiest case to consider. It includes 
materials such as the soft and medium hard metals. 
If the collision velocity is not too high, the damage 
marks produced on these materials by high-speed 
collision with liquid drops are empty spheroidal seg- 
ments. Surface and cross-sectional views of such 
pits in collision experiments performed elsewhere 2 
are shown in figures 1, 2, and 3. 

The damage pits produced in the planar surface of 
soft and medium hard metals when the leading sur- 
face of the solid runs into a liquid drop at high speed 
are almost exclusively due to the impact pressure 
produced. Although the radial flow of the liquid 
does make a contribution to the damage in the case 
of the soft metals (it drags metal up the walls of the 
damage pit and piles it up at the mouth of the crater), 
this contribution may be small enough in comparison 
with the damage that is caused by the pressure so 



that it can be neglected in first approximation. 
Neglecting all the dimensionless products that result 
from consideration of the radial flow of the liquid 
drop and from resilience of the solid target material, 
the approximate dimensionless damage equation for 
this class of materials is 

d'ld=F'[p'IP, c'/c, V/c, E'/(pV*)], (8) 

where 8' is the depth of the damage pits. 

2.4. Experimental Observations 

Graphs of measured values of pit depth plotted 
against relative collision velocity indicate that in 
the velocity range that has been investigated the 
pit-depth-versus-velocity curve is very close to a 
straight line regardless of whether the pits were 
caused by firing steel spheres against a stationary 
specimen plate or by firing a specimen plate at high 
velocity against relatively stationary liquid drops. 
It is to be expected that at a fixed relative collision 
velocity the depth of pit caused by a projectile that 
does not flow during the collision (hardened steel 
sphere) will be deeper than that caused by a projec- 
tile that does flow during the collision (liquid drop). 
In the first case all of the kinetic energy (neglecting 
that which is converted into heat) is delivered to the 
solid target; in the second case part of the kinetic 
energy is used to produce the flow of the projectile 
and only the remainder is delivered to the solid 
target. 

Experiments in which 7/32-in. steel spheres were 
fired against 1/8-in. -thick type-1100 aluminum 
plates have been carried out by Mr. Herschel L. 
Smith at the Naval Research Laboratory, Washing- 
ton, D.C Some of the data obtained in these ex- 
periments are presented in figure 4. One important 
fact is apparent from these data, namely, for the 
sphere diameter and for the thickness of specimen 
plate used, if the specimen plate has only peripheral 
support the pit-depth-versus-velocity curve has 
an intercept on the velocity axis, but if the specimen 
plate is rigidly backed by a heavy steel supporting 
plate the pit-depth-versus-velocity curve goes directly 
to the origin. 

Experiments in which target plates of different 
metals were fired into drops of mercury and drops 
of water at high speed have been carried out in an- 
other laboratory. 2 Some of the data obtained are 
plotted in figures 5 and 6. It can be seen that in 
every case the pit-depth-versus-velocity curve is 
essentially a straight line with an intercept on the 
velocity axis. From the bulge on the trailing face 
of the specimen in picture 6 of figure 2 it is apparent 
that the metal plates are supported during the 
collision with the liquid drops in such a way that 
the trailing face of the target plate is a free surface. 
An important fact apparent from the data plotted 
in figure 6 is that the slope of the straight line is 
a function of the drop size. 

On the basis of the evidence presented, if the speci- 
men has peripheral support only, the pit-depth- 
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I SURFACE VIEW 4 CROSS - SECTIONAL VIEW 

IMPINGEMENT VELOCITY 686 ft/sec 





2 SURFACE VIEW 5 CROSS -SECTIONAL VIEW 

IMPINGEMENT VELOCITY I250ft/sec 




3 SURFACE VIEW 6 CROSS— SECTIONAL VIEW 

IMPINGEMENT VELOCITY 1645 ft/sec 

Figure 1. Mercury-drop damage pits in lead. 
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I SURFACE VIEW 4 CROSS - SECT IONAL VIEW 

IMPINGEMENT VELOCITY 686 ft/sec 




2 SURFACE VIEW 5 CROSS - SECT ION AL VIEW 

IMPINGEMENT VELOCITY 1285 ft/sec 




3 SURFACE VIEW 6 CROSS -SECTIONAL VIEW 

IMPINGEMENT VELOCITY 2320ft/sec 

Figure 2. Mercury-drop damage pits in 1100 aluminum. 
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SURFACE VIEW 4 CROSS - SECTIONAL VI E W 

IMPINGEMENT VELOCITY 695 ft/sec 




2 SURFACE VIEW 5 CROSS - SECTIONAL VIEW 

IMPINGEMENT VELOCITY 1200 ft/sec 
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3 SURFACE VIEW 6 CROSS - SECTIONAL VIE W 

IMPINGEMENT VELOCITY 2445 ft/sec 

Figtjke 3. Mercury-drop damage pits in copper. 
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Figure 4. Best-fit curves for collisions of 0.55-cm steel spheres 
against plates of 1100 aluminum. 

O aluminum plate with edge support only; A aluminum plate rigidly clamped 
to heavy steel plate. 
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Figure 5. Best- fit curves for collisions of mercury drops 
against plates of three different metals. 

A 0.10-cm mercury drops against lead; □ 0.10-cm mercury drops against copper; 
O 0.15-cm mercury drops against iron. 

versus-velocity equation for materials of the response 
type being considered appears to be: 



h>=K x V-K 2i 

where K x and K 2 are constants. 5 



(9) 
It appears that 
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Figure 6. Best-fit curves for collisions of mercury drops of 
two sizes against plates of 1100 aluminum. 

O 0.10-cm mercury drops; A 0.20-cm mercury drops. 



this equation applies regardless of whether the pits 
are produced by the impingement of solid projectiles 
against solid targets or by firing solid target plates 
into liquid drops at very high velocities. However, 
because the amount of the kinetic energy delivered 
to the target depends upon whether or not the pro- 
jectile flows during the collision, the constants 
Ki and K 2 will be different for solid-target-versus- 
solid-projectile and for solid-target-versus-liquid- 
projectile collisions. 

2.5. Pit-Depth- Versus- Velocity Equation for Medium- 
Thin Metal Targets 

For targets that are metal plates, that have a 
thickness several times the diameter of the projectile, 
and that are mounted so that the reverse side of the 
plate is a free surface, a pit-depth- vcrsus-velocity 
equation can be developed by considering the move- 
ment of a cylindrical core of material through the 
target plate under the area of contact involved during 
the collision. If the target plate is fired against the 
projectile, this core of target material is slowed down 
with respect to the remainder of the target plate as 
a result of the collision and there is a relative motion 
between this core and the remainder of the target 
plate, which moves forward with respect to it. 



5 Partridge, VanFleet, and Whited [11] also reported a linear relationship be- 
tween pit depth and velocity; the linear relationship was established independ- 
ently by the author on the basis of the Convair data. 
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If the reverse side of the target plate is a free 
surface and the core can move freely with respect 
to the remainder of the plate, the kinetic energy 
transformed as a result of the collision is largely 
converted into the work of moving the core against 
the bonding that holds it to the remainder of the 
target plate. For these conditions there is an 
intercept on the velocity axis. In terms of this model, 
it would appear that the intercept velocity is related 
to the shear yield strength of the material of the 
target plate. In the range of velocities that are just 
in excess of the threshold velocity required to produce 
any permanent damage at all (intercept velocity), it 
appears that the core may only be displaced with 
respect to the remainder of the target plate. For 
this case plastic deformation occurs only at the 
cylindrical boundary of the core. At higher velocities 
plastic deformation may occur within the core itself. 

Consider that a small target plate moving at 
velocity V in a stationary coordinate system, which 
is located outside of the target plate, strikes a 
stationary liquid drop (see fig. 7). The result of the 
collision is that a wave of compression is initiated 
both in the solid material of the target plate and in 
the liquid of the drop. It is helpful to view the 
collision incident as it is seen by an observer located 
at the origin of the coordinate system in figure 7. 
This observer sees a target plate, all parts of which 
are moving forward uniformly at velocity V in the 
(+?/)-direction, approach a spherical drop, all parts 
of which are stationary. After the collision has 
occurred, the observer sees a zone in which the 
particles of the material of the target plate have 
taken on a velocity, aV, in the (— ^-direction where 
a is a constant; this zone is indicated schematically 
with a dotted line and the letter A in figure 7. 
Similarly, the observer sees a zone in the liquid drop 
in which the particles of the liquid have taken on a 
velocity, bV, in the (+y)-direction where b is 
a different constant; this zone is indicated schemati- 
cally with a dotted line and the letter B in figure 7. 
Zone A spreads at the characteristic sound velocity 
of the target material, c', through the thickness, d' , of 
the target plate and zone B spreads at the charac- 
teristic sound velocity, c, of the liquid of which the 
drop is composed, through the diameter, d, of 
the liquid drop. 

To the observer who looks first at a point A' that 
is in the target plate but outside the spreading 
boundary of zone A and who then looks at the mate- 
rial within zone A, it appears that the particles at 
A' have a velocity V and that the particles in zone A 
have a velocity (l-a)V. After a time t that is just 
long enough for zone A to complete one trip through 
the thickness of the target plate, the leading surface 
of the target plate at point A' in figure 7 has ad- 
vanced a distance Vt, but in the same time t the 
area of the leading surface of the target plate that 
constitutes the upper boundary of zone A in figure 
7, has advanced only by {\-a)Vt. If the average 
negative velocity of the core is taken to be aV/2, 
the displacement that exists between the leading 




Figure 7. Stationary coordinate system for a liquid-drop- 
versus-solid-specimen collision. 

surface of the target plate at point A' and the upper 
boundary of zone A after time t is aVt/2. 

If the observer watches zone B as it just spreads 
to the opposite side of the liquid drop, he sees that 
this zone reflects from the free air-liquid interface 
with change of sign, that is, he sees a tension wave 
begin to move toward him in the (—y) -direction 
of figure 7. Within the tension wave the particle 
velocity is 2bV in the (+2/) -direction. The tension 
wave moves toward the observer at the characteristic 
speed of sound, c, of the liquid of which the drop is 
composed. In the zone of the liquid of the drop 
that has been traversed by the tension wave the 
particles of the liquid are unstressed and are moving 
at velocity 2bV in the (+2/) -direction. At the 
instant that the tension wave returns to the collision 
surface, where radial flow of the liquid of the drop 
is occurring, all of the liquid of the drop is unstressed 
and moving at the velocity 2bV in the (-\-y)~ 
direction. 

It can be assumed that return to the collision 
surface of the tensile wave in the liquid provides a 
cutoff for the collision. At this instant all of the 
liquid of the drop is unstressed. The collision sur- 
face will appear to be a free surface to the com- 
pressed particles of the material of the target plate 
in zone A. A wave of relief (tension) will then be 
initiated in the material of the target plate and will 
begin to trail zone A which continues to spread 
through the thickness, d' , of the target plate. 
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Let the arbitral time t referred to previously 
now be the time that is required for zone B to make 
one round trip through the liquid drop. Then 
t=2d/c. The distance to which zone A has spread 
in this time is c't and this determines the length of 
the compressed zone A wdiich continues to move 
through the target plate. When zone A reaches the 
free trailing surface of the target plate, it also 
reflects as a tension wave. This tension wave moves 
in the (+y) -direction of figure 7. In that part of 
zone A that has been traversed by the returning 
tension wave the particles of the material of the 
target plate are under zero stress and have velocity 
(1 — 2a) V in the (+y) -direction. When the re- 
flected tension wave in the target plate reaches the 
leading surface of the plate it undergoes partial 
reflection back into the material of the target plate 
and partial transmission into the liquid of the drop 
which is running off radially around the central point 
of the collision. 

a. Slope of the Pit-Depth- Versus- Velocity Curve 

The compressional wave spreads slightly as it 
moves through the thickness of the target plate; 
consequently, the core of material through the target 
plate that is slowed down as a result of the collision 
is in reality somewhat conical in shape. For sim- 
plicity, the true situation may be idealized in two 
ways. First, the core may be regarded as a true 
cylinder which is free to move or slide in the direc- 
tion of the collision blow but which is restrained 
laterally. A similar cylinder exists in the liquid of 
which the drop is composed. Secondly, the compres- 
sional waves that move through the cylinder in the 
target plate and through the cylinder in the drop 
may be regarded as plane waves. 

For the case of plane waves, a=z/(z-\-z f ) and 
b = z'l(z+z') where z=cp and z' = c'p are the 
acoustic impedances of the material of the drop and 
of the material of the target plate, respectively. 
The relations are found by equating the plane-wave 
stress in the material of the drop, o-, to the plane-wave 
stress in the material of the target plate, g\ that is, 



-cpv=c f pv f - 



(10) 



where v, v r are the particle velocities in the zones 
traversed by the compressional waves initiated in 
the drop and in the target plate, respectively, and 
by imposing the condition that the impacting surfaces 
remain in contact, namely, 



v+v'=V. 
From eq (10) and (11), 

f>'=c P T7(cp+cV) 

= zV/(z+z') 
=aV 



(11) 



niid 



v=z'VI(z+z') 
=bV. 



(12) 



This derivation is the same as that for the collision 
of two free rods except that, for the case that the 
rods are cylinders that are free to move in the colli- 
sion direction but that are restrained laterally, the 
sound speed that must be used is the speed of sound 
in an infinite medium. This can be shown as 
follows. 6 According to Hooke's law: 



Eti =ai —v(2a r ) = ai —2v<T r , 



(13) 



where E is Young's modulus, e t is the longitudinal 
strain, <n is the longitudinal stress, v is Poisson's 
ratio, and o T is the radial stress. Also, 

Ee r =a r -p(a l + ( r r )=0 y (14) 

where e r is the radial strain. From eq (14), 

<r r =wr,/(l-p). (15) 

By substituting the expression for a r given by eq 
(15) into eq (13) it is found that 



{l-v)Eu 

(Ji= (i-2v){\+vy 



(16) 



and because, for plane waves, c=[ai/pei] 1/2 

This is the sound speed for an infinite medium and 
it is the sound speed that must be used both for c 
and for c' in the expressions for a and for b. 

The depth of pit produced at any relative impinge- 
ment velocity is proportional to the negative velocity 
produced in the cylindrical core in the target plate 
as a result of the collision and to the time that this 
negative velocity exists. If it is assumed that the 
elastic wave that is induced in the cylindrical core 
as a result of the collision makes a round trip through 
the core in the time interval 2d/c, then, at the end 
of this time interval, the cylindrical core is moving 
as a rigid body and the negative velocity that it has 
acquired with respect to the remainder of the target 
plate is given approximately by 2zV/(z+z'). During 
this time interval the average negative velocity given 
to the material of the cylindrical core was zV/(z+z'). 
The slope K x of eq (9) is then k x (2d\c) [z/(z+z f )] 
where ki is a new constant. 

The experimental data for a number of target 
metals of widely different properties fired against 
drops of mercury, and for three target metals fired 
both against drops of mercury and drops of water, 
require that ki = 3.Q. Because this value of the 
constant Id brings the calculated curve into good 
agreement with the experimental data for a number 
of target metals and for two drop liquids, it would 
appear that k x does not involve properties either of 

e The author is indebted to Dr. John M. Frankland of NBS Mechanics Section 
for the proof in elasticity theory that it is the sound speed in infinite medium that 
is required. 
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the material of the target plate or of the liquid of the 
drop. Therefore, eq (9) is 

Substitution of z=cp, z' = c'p', and division by d puts 
eq (18) in the dimensionless form 



5' 



7.2 



1 + 



am 



©- 



■K 2 /d, 



(19) 



where K 2 /d is a dimensionless intercept. Four of the 
dimensionless quotients predicted by eq (8) are found 
in eq (19). 

b. Intercept on the Velocity Axis 

The process previously postulated, namely, that 
the core of target material below the collision area 
is slowed down with respect to the remainder of the 
target plate so that the remainder of the target plate 
moves forward with respect to it, certainly occurs 
at all impingement velocities. However, no ob- 
servable damage pit is produced at all below a 
critical velocity, and, if the liquid of the drop is 
kept constant, this critical velocity is different for 
each of the metals for which data are available. 
(See fig. 5.) The intercept velocity is also different 
on the same metal if the liquid of the drop is changed. 

The fact that no permanent damage pit is formed 
below a given velocity which is characteristic for 
each metal, indicates that for velocities lower than 
this critical velocity the relative motion between 
the core and the remainder of the target plate is 
completely elastic and that no permanent shear 
deformation occurs in the material around the 
cylindrical wall of the core. The intercept of the 
pit-depth-versus-velocity curve appears to be a 
function of the shear yield strength of the material, 
that is, of the remaining dimensionless quotient 
E'/(pV 2 ) of eq (8) taking E' ', the energy per unit 
volume absorbed by the target material without 
fracture or plastic yield, to be the shear yield 
strength. 

It has been pointed out that the experimental 
pit-depth-versus-velocity curve is a straight-line 
function of the velocity. Therefore, the dimension- 
less quotient E' l(pV 2 ) cannot be used for the inter- 
cept without eliminating the factor 1/V 2 . This can 
be accomplished by multiplying the dimensionless 
quotient E'/(pV 2 ) by the square of the dimensionless 
quotient V/c of eq (8) to obtain E f j(pc 2 ). If this 
expression is substituted for the dimensionless inter- 
cept K 2 /d in eq (19) and if the intercept condition 
that 8 f /d=0 is imposed, the expression that is found 
for the intercept velocity is not able to account for 
the observed ratio of experimental intercepts. It was 
found that the observed intercept velocities can be 
accounted for if the dimensionless quotient E f /(pc 2 ) 
is multiplied by the dimensionless quotients 
(p/p') 1/2 > z/z't c / c> ', an( i by a dimensionless numerical 



constant, k 2 , having the value 136.8. Substitution 
of the resulting dimensionless quotient, 136.8 
E'z/(p 1/2 p n/2 cc'z'), for the dimensionless intercept, 
K 2 jd, in eq (19) and use of the acoustic impedance 
z for cp produces the dimensionless pit-depth equa- 
tion 

8' 7.22_ /V\ 136.8ff'z 

■ ' ^ C J ~W~'l/2 M f„f' 



d z+z' 



pWp' 1 ' 2 ' 



cc z 



If the intercept condition that 8'jd=0 is imposed, 
the intercept velocity, V iy is found to be given by 



v<= 



WE' (z+z') 



(sc'z'y* 



(21) 



It is noteworthy that one would expect V t to be 
given by E r (z+z')/zz' multiplied by some constant, 
because the plane wave stress is given by zz'Vf 
(z-\-z'). It was found, however, that this expected 
expression would not give the observed intercept 
velocities for both mercury drops and waterdrops if 
the same numerical constant was used for both drop 
liquids. Equation (21) appears at present to be the 
best expression for the intercept velocity. When 
more data are obtained and the problem is studied 
further, it may be found necessary to change it. 

The numerical constants in eqs (20) and (21) have 
been chosen to give best fit to the pit-depth-versus- 
velocity data. The values used for the physical 
constants of the different materials are given in 
table 2. In choosing the constants, the speed of 
sound in infinite medium was used for c' and the 
dynamic compressive yield strength was used for E'. 
(In liquids there is only one speed of sound that can 
be used for c.) Dynamic rather than static strength 
values must be used for E' because the loading time 
is of the order of several microseconds. The strength 
that should be used for E' is the dynamic yield 
strength in shear. Unfortunately, the amount of 
work that has been done in determining the dynamic 
yield strengths of materials is very limited. Whiffin 
[7] has determined the dynamic yield strength in 
compression for a number of materials. In choosing 
the numerical constant for eq (21), the dynamic 
yield strength in compression was used for E r in- 
stead of the dynamic yield strength in shear. Some 
justification for this substitution may be found in 
the Von Mises strain energy theory according to 
which the tensile elastic limit, <r y , is V3 r y where r y is 
the yield strength in shear. 

2.6. Test of the Equation 

For pit-depth calculations, eq (20) may be put 
in the form 



°-m?) iv - v > ] ' 



(22) 



where V t is the intercept velocity for the particular 
liquid-drop-solid-target combination being used and 
where V t is given by eq (21). Pit-depth-versus- 
velocity data obtained in another laboratory 2 were 
used to test the equations. In most cases the 
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Table 


2. Physical Constants 






Metal 


Density, p 


Sound speed, c 
(infinite medium) 


Acoustic imped- 
ance, 2 


Static tensile 

yield 
strength, Y 


Dynamic compres- 
sive yield 
strength, E' 


1100-0 aluminum. 


g/cm* 
* 2. 713 
a 2. 768 

b8.92 
b 11. 3437 


cmjsec 
d 6. 318X105 
d 6. 370X105 

d 4. 691X105 
d 2. 277X105 
d 2. 128X105 
d 5. 786X105 
e 5. 850X105 

e 4. 170X105 

e 1.451X105 
el. 497X105 


g /sec- cm 2 
1.714X106 
1.763X106 

4. 184X106 
2.583X106 
2.414X106 
4. 547X106 
4.598X106 
2.977X106 
1. 966X106 
0. 1493X106 


psi 
/ ' 2, 625 
I o 5, 000 
* 12, 625 

' 3, 975 


d/cm* 
k 7. 239X108 
fa 1. 183X109 
« 2. 350X109 

fa 2. 394X109 
1*4.324X108 


2024-0 aluminum _ _ 


(Nipper, annealed electrolytic tough 
pitch _. 


Lead, chemical ___ 


Lead, pig __. ._ _ . 




Steel, cold rolled 


e 7. 859 
b7.86 
b7. 14 
b 13. 546 
b 0. 99707 




* 8. 274X109 
h 7. 660X109 


Iron 




Zinc 




Mercury 






Water, 25° C. 













» Data from Aluminum Co. of America. 

b Data from Handbook of Chemistry and Physics. 

o Data from Metals Handbook, 1948. 

d Measured in NBS Sound Section by Carroll Tschiegg. 

e Data from L. Bergmann [12]. 

physical properties of the metals used in the experi- 
ments were not determined. The physical con- 
stants used for these materials in eqs (21) and (22) 
are given in table 2. 



a. Test of the Equation With Mercury Drops and Pure Metal 
Targets 

Pit-dept h-\ nsus-velocity curves calculated by 
eq (22) for collisions of 0.10- and 0.15-cm mercury 
drops against copper are shown in figure 8. Experi- 
mental points for the two drop sizes are indicated in 
this figure with circles and with triangles, respec- 
tively. The material used in these experiments 2 was 
described as pure copper. Physical constants for 
electrolytic copper were used in the calculations. 
The intercept velocity, V iy was calculated by use 
of eq (21). It can be seen that both the slope and 
intercept of the calculated curves are in good agree- 
ment with the experimental data. 

Pit-depth-versus-velocity curves calculated by 
use of eq (22) for collisions of 0.10-, 0.15-, and 0.285- 
cm mercury drops against 1100-0 aluminum are 
shown in figure 9. Experimental points for the 
three drop sizes that were used are indicated in this 
figure with circles, triangles, and squares, respec- 
tively. The intercept velocity, V i} was calculated 
by use of eq (21). The value of the dynamic com- 
pressive yield strength that was used for E f in eq 
(21) was calculated by means of a formula given 
by Whifnn [7] for the Duralumins. The static 
yield strength of the 1100-0 aluminum used was not 
known. 2 For the purpose of the calculation of E' it 
was taken to be 5,000 psi (table 2). The speed of 
sound in infinite medium, determined for a piece 
of 1100-0 aluminum of lower yield strength by 
Tschiegg (table 2) was used for c' and in computing 
z f . There is more scatter in the experimental data 
for the 0.15-cm drop size than for the other two 
drop sizes. Scatter could be caused by variation 
of the drop size from the reported value; it could 
also be caused by use of target plates cut from sheet 
material of different yield strength. The data for 
both the 0.15-cm and 0.285-cm drop sizes would fit 
the curves better if the intercept velocity were 
1.1X10 4 cm/sec rather than 1.26X10 4 cm/sec, which 



' Measured in NBS Engineering Mechanics Section by Lafayette 
K. Irwin, 
g Data of A. C. Whiffin, communicated by letter, 
fa Data of A. C. Whifhn [7]. 
» Extrapolated from graph, figure 7e, given by Krafft and Sullivan [9]. 

was the value found by use of eq (21). The experi- 
menters report [8] that the target plates were cut 
from 1100-aluminum sheets having thicknesses in 
the range of 0.185 to 0.210 in.; hence, they were not 
all cut from the same sheet stock and, consequently, 
some variation in yield strength, which would affect 
the intercept velocity directly, could have been 
present. In general, however, there is good agree- 
ment between the experimental data and the 
calculated curves. 




3 4 5 6 7 

IMPINGEMENT VELOCITY , cm /sec 



Figuee 8. Calculated curves for collisions of mercury drojs of 
two sizes against plates of copper. 

O observed depth for 0.10-cm drop; A, calculated curve for 0.10-cm drop; 
A observed depth for 0.15-cm drop; B, calculated curve for 0.15-cm drop. 
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Figure 9. Calculated curves for collisions of mercury drops of 
three sizes against plates of 1100 aluminum. 

O observed depth for 0.10-cm drop; A, calculated curve for 0.10-cm drop; A 
observed depth for 0.15-cm drop; B, calculated curve for 0.15-cm drop; El observed 
depth for 0.285-cm drop; C, calculated curve for 0.285-cm drop. 



.24 




1 




1 




1 
A 




' B 




.22 


_ 


















.20 


- 
















- 


. 1 8 


- 








A 








" 


. 1 6 








L 










- 


£.14 


















- 


o 




















I 




















CL .12 






A 












- 


Ul 

o 




















H 


















o 


CL . 1 


















A 


.08 


- 












o 




- 


.06 


- 








Oy 


O ^ 






- 


.04 


~ 

















" 


.02 


- 


/ O/ 


Orf^ 












- 







L^ i 




1 




1 




l 





2 3 

IMPINGEMENT VELOCITY ,cm/sec 



5XI0 4 



Figure 10. Calculated curves for collisions of mercury drops 
of two sizes against lead target plates. 

O observed depth for 0.10-cm drop; A, calculated curve for 0.10-cm drop; A ob- 
served depth for 0.285-cm drop; B, calculated curve for 0.285-cm drop. 



The pit-depth-versus-velocity curves calculated 
by use of eq (22) for collisions of 0.10-cm and for 
0.285-cm mercury drops against lead are shown in 
figure 10. The experimental points for the two drop 
sizes are indicated with circles and with triangles, 
respectively. The intercept velocity, V iy was cal- 
culated by use of eq (21). The value of the dynamic 
compressive yield strength used for E' in eq (21) 
was that given for 99.97 percent lead by Whiffin 
(table 2). The value of the speed of sound in infinite 
medium for pig lead measured by Tschiegg (table 2) 
was used for c' and in computing z' . The calculated 
curves are a good fit to the data at the lower veloci- 
ties; at the higher velocities the depth of the pits is 
greater than that predicted by eq (22). The pits 
produced in lead by mercury drops are characterized 
by a heavier lip of metal around the mouth of the 
crater than is observed on pits in 1 100-0 aluminum or 
on pits in copper. See figures 1,2, and 3. It ap- 
pears that for a metal that is as soft as lead the flow 
of a liquid-drop projectile at high collision velocities 
drags a notable amount of metal from the bottom 
of the pit and piles it up at the mouth of the crater. 
This extra mode of pit formation was neglected in 
the development of eq (22) and, therefore, this 
equation does not fully account for the depth of 
pits formed at very high velocities in metals as soft 
as lead. 



b. Test of the Equation With Mercury Drops and Alloy Metal 
Targets 

Experimental pit-depth-versus-velocity data were 
obtained for collisions of 0.15-cm mercury drops 
with steel targets in another laboratory. 2 It was 
reported that the target material was originally 
%H cold-rolled steel. After annealing above the 
alpha temperature it developed a case hardening. 
The case hardening was chipped off and the target 
plates were polished. It was reported that the 
resulting material had an average Rockwell E hard- 
ness of 90.96 and Brinell hardness of 90.98, which 
they stated corresponds to an ultimate tensile 
strength of 47,200 psi. 

In order to calculate the intercept velocity by 
means of eq (21) it is necessary to know the dynamic 
compressive yield strength of the target metal. 
The manganese and carbon contents of steels strongly 
affect their dynamic yield strength; the manganese 
content is particularly important if it is low, for 
example, in the range from zero to 0.50 percent. A 
chemical analysis of the manganese and carbon 
content of one of the steel target plates was made 
in the NBS Analytical Chemistry Section. It re- 
vealed that the steel contained 0.70 percent man- 
ganese and 0.21 percent carbon. The value of the 
dynamic compressive yield strength used for E' in 
eq (21) was found by extrapolating a curve given by 
Krafft and Sullivan [9] for a steel containing 0.98 
percent manganese and 0.22 percent carbon that was 
loaded in compression. The speed of sound in 
infinite medium determined for cold-rolled steel by 
Tschiegg (table 2) was used for c f and in com- 
puting z' . 
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The pit -depth-versus- velocity curve calculated by 
use of eqs (21) and (22) for collisions of 0.15-cm 
mercury drops against steel targets is shown in 
figure 11 along with the experimental points. It 
would appear that the intercept velocity is about 
correct. However, a straight line that would fit 
the experimental data would have a slope consider- 
ably less than that of the calculated curve. It 
would appear from this evidence that eq (22) does 
not predict pit depth with exactness for pits formed 
in targets of alloy steel. 




3 4 5 6 7 

IMPINGEMENT VELOCITY cm/sec 



Ficure 11. Calculated curve for collisions of 0.15-cm mercury 
drops against steel. 

O observed depth for 0.15-cm drop. 



The pit-depth-versus-velocity curves calculated 
by use of eqs (21) and (22) for collisions of 0.10-cm 
and 0.20-cm mercury drops against 2024-0 aluminum 
targets are shown in figure 12 along with experimental 
points obtained elsewhere. 2 The experimental points 
for the two drop sizes are indicated with circles and 
triangles, respectively. The targets were prepared 
at the NBS. The average static yield strength 
(0.2% offset) of the metal used was found to be 12,625 
psi from measurements made by Irwin (table 2). 
The dynamic compressive yield strength used for 
E' in eq (21) for this material was calculated from a 
ratio of the static to the dynamic yield strength of 
2.7. This ratio was estimated by W'hiffm [10] for the 
2024-0 aluminum that was used for the targets. 
The speed of sound in infinite medium determined 
for 2024-0 aluminum by Tschiegg (table 2) was used 
for c' and in computing z* '. The calculated curves 
differ from a best-fit line drawn through the experi- 
mental points in intercept but not in slope. If the 
intercept velocity had been 2xl0 4 cm/sec rather than 
2.42 lxl 4 cm/sec there would have been relatively 
good fit between the calculated curves and the 
experimental points. 

In summarizing, it can be said that eqs (21) and 
(22) have not given calculated curves that show as 




4 5 6 7 

velocity, cm/sec 

Figure 12. Calculated curves for collisions of mercury drops 
of two sizes against target plates of 202J+-0 aluminum. 

O observed depth for 0.10-cm drop. 

A, calculated curve for 0.10-cm drop. 
A observed depth for 0.20-cm drop. 

B, calculated curve for 0.20-cm drop. 

good agreement with experimentally determined 
points for collisions of target plates of alloy metals 
with mercury drops as that obtained for collisions 
of plates of pure metals with mercury drops. 

c. Tests of the Equation With Waterdrops 

Test results with waterdrops of constant size are 
available for some of the metals used in the experi- 
ments with mercury drops. 

The pit-depth-versus-velocity curve calculated by 
use of eqs (21) and (22) for collisions of 0.2-cm 
waterdrops against targets of annealed electrolytic 
tough pitch copper is shown in figure 13 along with 
experimental points obtained elsewhere. 2 The tar- 
gets were prepared at the NBS. The static yield 
strength (0.2% offset) of the metal used was deter- 
mined by Irwin; the average value was 3,975 psi. 
The dynamic compressive yield strength used for 
E' in eq (21) was thai riven by Whiffin (table 2) 
for electrolytic copper. The speed of sound in in- 
finite medium determined for electrolytic tough 
pitch copper by Tschiegg (table 2) was used for c r 
and in computing z f . There is quite a bit of scatter 
in the data plotted in figure 13. However, both 
the slope and intercept of the calculated line are in 
relatively good agreement with these data. 

The pit-depth-versus-velocity curve calculated by 
use of eqs (21) and (22) for collisions of 0.2-cm 
waterdrops against targets of 2024-0 aluminum is 
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shown in figure 14 with experimental points obtained 
elsewhere. 2 The targets were prepared at the NBS. 
They were made of the same 2024 aluminum as the 
targets that were used for collisions with mercury- 
drops (see sec. 2.6b and fig. 12), and the same 
values of E', c', and z f were used in calculating the 
pit-depth-versus-velocity curve. It can be seen 
from figure 14 that the calculated line fits the experi- 
mental points fairly well. 
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Figure 13. Calculated curve for collisions of 0.20-cm water- 
drops against target plates of electrolytic tough pitch copper. 

G observed depth for 0.20-cm drop. 

The pit-depth-versus-velocity curve calculated by 
use of eqs (21) and (22) for collisions of 0.2-cm 
waterdrops against leac' targets is shown in figure 15. 
The experimental points shown in the graph were 
obtained elsewhere. 2 The kind of lead that was 
used for the targets is not known. The dynamic 
yield strength used for E' w^as that given bv Whiffin 
(table 2) for 99.97 percent pure lead. The speed of 
sound in infinite medium for pig lead determined by 
Tschiegg was used for c f and in computing z\ 

The amount of agreement that has been found 
between the calculated curves and experimental 
points for collisions of both mercury drops and water- 
drops with targets of copper, 2024-0 aluminum, and 
lead may be favorable evidence for the possible 
usefulness of eq (22) in predicting corresponding 
velocities-for-equal-pit-depth for collisions of metal 
target plates with waterdrops and with mercury 
drops. 

2.7. Corresponding Velocities-for-Equal-Pit-Depth 

From the standpoint of the problem of high-speed 
rain erosion, the practical value of an equation giving 
the depth of damage pits as a function of impinge- 
ment velocity is its use to extrapolate from the 
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Figure 14. Calculated curve for collisions of 0.20-cm water- 
drops against target plates of 202^-0 aluminum. 

O observed depth for 0.20-cm drop. 
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Figure 15. Calculated curve for collisions of 0.20-cm 
waterdrops against target plates of lead. 

O observed depth for 0.20-cm drop. 



known depth of pit produced by collision of a solid 
target with a drop of high-density liquid at a rela- 
tively low impingement velocity to the depth of pit 
that would be produced on the same target material 
by collision with a drop of low-density liquid, such 
as water, at a very much higher impingement ve- 
locity. The purpose of carrying out this extrapola- 
tion is to bypass the necessity of firing target mate- 
rials at the extremely high velocities for which test 
results are desired for collisions of solids with water- 
drops. 

If drops of a liquid A having diameter d A collide 
with a solid material having acoustic impedance z' 
and dynamic yield strength E' ', eq (22) predicts that 
the pit depth will be 
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«*'-^SV r '- ViA 



(23) 



where V iA is the intercept velocity calculated with 
use of eq (21) and with use of the density and acoustic 
impedance of liquid A for p and z, respectively. The 
depth of pits resulting from collisions of drops of 
liquid B with the same solid material will be 



*i= 



lM B z L 



C B {z B + Z f ) 



[V B -V iB ], 



(24) 



where V iB is the intercept velocity calculated with 
use of eq (21) and with use of the density and acoustic 
impedance of liquid B for p and z, respectively. 

If the depth of pits produced by collision of the 
solid target with drops of liquid A is to be the same 
as that produced by collision of the solid target with 
drops of liquid B, then, bv equating 8 A and b B given 
by eqs (23) and (24), 



V A 



_ d B z B c A (z A +z') 
~d A z A c B (z B +z') 



Vb 



d B z B c A (z A +z')p9E' (z B +z') l 



d A z A c B (z B +z')L ( PB c'z'*y>* 



19E'(z A + z') 
' (pac'z'*) 



,'3\l/2 



(25) 



All quantities are expressed in cgs units. The dy- 
namic compressive yield strength must be used for E f 
and the sound speed in infinite medium must be 
used for c' and in computing z f in making calcula- 
tions with this equation. 

2.8. Applicability of the Damage Equation 

It was remarked in section 2.2 that the extent and 
type of damage produced on solids as a result of col- 
lision with liquid drops depends strongly on the 
characteristic properties of the solid, and that if the 
damage is to be described by an equation, it will be 
necessary to develop separate equations for the 
damage produced on solid materials that have widely 
different properties. It was pointed out that the 
simplest case to consider is that of materials such as 
the soft and medium hard metals that undergo 
permanent plastic flow without fracture as a result 
of collision with liquid drops. Equation (22) has 
been developed to describe the depth of pits that will 
be produced in materials of this kind as a result of 
such collisions. It cannot validly be applied beyond 
the limits of the model on which it was constructed. 

Equation (22) is restricted to damage associated 
with the impact pressure that results when a solid 
target collides with a liquid drop at high speed. 
The dimensionless Reynolds number, Weber num- 
ber, and surface tension quotient, which are associ- 
ated with the radial flow r of the drop, were deleted 
from eq (7) to produce eq (8), and eq (22) contains 
only the dimensionless quotients that appear in eq 
(8). In view of this restriction on eq (22) it can be 



expected that there will be limits on its application 
even to the soft and medium hard metals. 

In discussing the agreement of the calculated 
curve with the experimental points for lead target 
plates it was pointed out that for metals as soft as 
lead at very high collision velocities the radial flow 
of a mercury drop drags metal up the walls of the 
damage pit and piles it up at the mouth of the crater. 
This mode of forming, or of deepening, a pit was 
not considered in the development of eq (22), and 
eq (22) will not adequately describe pits that Avere 
produced either wholly or appreciably by this mech- 
anism. It is noteworthy that where tins additional 
pit-forming mechanism does not operate (in the low- 
velocity range), eq (22) does adequately describe 
the pits that are formed in lead. 

Two other points should be mentioned with regard 
to the valid use of eq (22) for de termini nir the depth 
of pits that will be produced in collisions between 
target plates of the medium hard metals (to which it 
applies) and liquid drops. One is the mode of 
mounting the plate of solid material with which 
the liquid drops collide; the other is the thickness 
of this plate. The reverse side of the plate must 
be maintained as a free surface for otherwise the 
core of solid material through the metal plate under 
the contact or collision area cannot move freely with 
respect to the remainder of the plate, which is an 
essential of the model on which eq (22) is based. 
The bulge on the reverse side of the 1100-aluminum 
plate shown in picture 6 of figure 2 appears to indi- 
cate that this condition was sufficiently realized in 
the mounting of the metal plates that were used for 
the mercury drop experiments. 2 If this condition is 
not realized, eq (22) cannot be applied to the pit 
depths that are obtained. The thickness of the plate 
should not be less than 1.5 to 2.0 times nor greater 
than 4 to 5 times the diameter of the impinging 
liquid drops. If the plate is too thin, tin 4 model on 
which eq (22) is based will break down because the 
plate will bend as a unit under the collision blow. 
If the plate is too thick, the spreading of the com- 
pressional wave on passing through it may cease to 
be negligible. 

Equation (22) was constructed on the assumption 
of permanent plastic flow of the material of the 
target plate. The dimensionless quotient e'/E' of 
eq (7), which represents the resilience of the plate 
material, has been neglected. Consequently, eq (22) 
cannot validly be applied to determine the depth of 
pits that will form in highly resilient materials such 
as the polymers and rubbers when they collide with 
liquid drops. In the case of such materials the core 
of target material depressed as a result of the 
collision tends to spring back into its original posi- 
tion; the permanent damage mark that remains 
after the collision is more nearly a circular cut than 
a pit. 

Equation (22) has been constructed using a very 
simple model. To determine how the Reynolds 
number and the Weber number should be introduced 
into eq (22), pit-depth-versus-velocity data should 
be obtained on a single metal for the condition that 
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the viscosity and surface tension of the arop liquid 
are gradually changed (such as by the use of glyeerol- 
water solutions). To determine how the sound 
speed, c, of the liquid should appear in the expression 
for the intercept velocity, data should be obtained 
using the same metals that have been used already 
and a liquid that has a sound speed different from 
those of water and mercury, which are nearly the 
same. Chloroform has a sound speed about 68 per- 
cent that of water. The sound speed of water can 
be increased by adding sodium chloride. The speed 
of sound in a 20 percent sodium chloride solution is 
about 13 percent greater than that in fresh water and 
it is indicated that increase in the salt concentration 
will further increase the sound speed. 

Finally, more pit-depth-versus-velocity data 
should be obtained for collision of both mercury 
drops and waterdrops against target metals whose 
dynamic compressive yield strengths are known. 
Whiffm [7] has determined the dynamic compressive 
yield strength of standard silver (7.5% Cu, 92.5% 
Ag), electrolytic copper, 99.97 percent lead ,and 
Armco iron (0.016% C, 0.006% Si, 0.017% S, 
0.003% P, and 0.030% Mn). In the process of 
verifying eq (22) it is important to use target ma- 
terials whose dynamic compressive yield strengths 
are known, and the same target material should be 
used throughout where drops of different liquids are 
used. If eq (22) is fully verified, it can be used to 
determine the dynamic compressive yield strength 
of other metals. 

3. Collisions Between Metal Plates and 
Flowing Metal Spheres 

Partridge, VanFleet, and Whited [11] fired spheres 
of zinc, tin, copper, lead, aluminum, and iron against 
targets of the same material at collision veloc- 
ities up to 24X10 4 cm/sec (7,900 ft/sec) and main- 
tained conditions such that the spherical pellet 
lost no mass before striking the target. They found 
that the penetration varied linearly with velocity for 
the materials used and the velocity range investi- 
gated. They reported that the pellets flowed during 
the collisions. 

In the light of the fact that the pellets flowed 
during the collisions, it is reasonable to suppose that 
the pit-depth-versus-velocity equation that was 
developed for collisions of metal targets with liquid 
drops, eqs (21) and (22), should apply to this case 
without modification. When the material of the 
target is the same as the material of the projectile, 
eqs (21) and (22) simplify to 



and 



S' = 3.6d(y-W 



V t =38 E'/z 



(26) 



(27) 



Some of the pit-depth-versus-velocity data [11] are 
plotted in figures 16, 17, 18, 19, and 20 along with 



curves calculated by use of eqs (26) and (27). The 
physical properties of the metals used by Partridge 
et al. [11] were not determined. The physical 
constants used for these materials in eqs (26) and 
(27) are given in table 2. 

3.1. Collisions of Iron Spheres Against Iron Targets 

It is assumed that soft, relatively pure iron was 
used in these experiments [11] because the spheri- 
cal pellets fired were made by placing fragments 
of the metal in the hemispherical cavities of a 
case-hardened steel tool and pressing the two 
sections together. The diameter of the pellets for 
which the experimental pit-depth-versus-velocit\* 
data are given in figure 16 was 0.483 cm. The 
calculated curve is shown in figure 16 along with 
the experimental points which are indicated with 
circles. The observed pit depths for iron spheres 
colliding with iron targets are in good agreement 
with the calculated curve up to a collision velocity 
of 24X10 4 cm/sec (7,900 ft/sec). This velocity is, 
however, less than half the speed of sound in infinite 
medium for iron (58.5 X10 4 cm/sec). 
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Calculated curve for collisions of 0.483-cm iron 
spheres against iron targets. 



FlGUKE 16. 

O observed depth for 0.483-em sphere, 



3.2. Collisions of Aluminum Spheres With Aluminum 
Targets 

Values for physical constants of 1100-0 aluminum 
were used in eqs (26) and (27). The value of the 
dynamic compressive yield strength used for E' 
in eq (27) was computed from a ratio of the dynamic 
to the static vield strength of 4; this is a rough 
ratio estimated by Whiffin [10] for an 1100-0 
aluminum having an average static vield strength 
of 2,625 psi. 

Very few experimental pit depths were given by 
Partridge et al. [11] for collisions of 0.483-cm alumi- 
num spheres with aluminum targets. The data 
available contain quite a bit of scatter. They are 
plotted in figure 17 along with the curve calculated 
using eqs (26) and (27). There is reasonably good 
agreement with the calculated curve. The highest 
collision velocity for which a pit depth was reported 
is approximately one-third of the speed of sound in 
aluminum. 
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Figure 17. Calculated curve for collisions of 0.483-cm 
aluminum spheres against aluminum targets. 

O observed depth for 0.483-cm sphere. 

3.3. Collisions of Lead Spheres Against Lead Targets 

Whiffin [7] determined the dynamic compressive 
yield strength of lead that was 99.97 per cent pure 
and this value of the dynamic compressive yield 
strength was used for E' in eq (27). The sound 
speed in infinite medium was determined by Tschiegg 
(table 2) both for chemical lead and for pig lead. 
The theoretical curve computed using the sound 
speed for pig lead is shown in figure 18 along with 
the pit-depth-versus-velocity data obtained by 
Partridge et al. [15] for 0.483-cm lead spheres 
impinging against lead targets. 
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Figure 18. Calculated curve for collisions of z /\&-in. lead 
spheres against lead targets. 

O observed depth for M6-in. sphere. 



The experimental pit depths found at the high 
velocities used are all greater than eq (26) would 
predict. No data were obtained in the low velocity 
range where the calculated curve was found to fit 
the data for collisions of lead target plates with 
mercury drops. See section 2.6a, where the pos- 
sibility that two mechanisms contribute to the 
formation of pits in lead at high collision velocities 
is discussed. 



The speed of sound in infinite medium for pig lead 
is marked with a vertical dotted line in figure 18. It 
would appear from the empirical data at the highest 
velocities used that the experimental pit-depth- 
versus-velocity curve for lead is flattening at veloci- 
ties above the speed of sound in lead. 

It appears from the data that have been presented 
in the preceding sections that, at collision velocities 
for which an originally solid target remains solid dur- 
ing the collision, pit depth is a straight-line function 
of collision velocity if the rear face of the target plate 
is a free surface. This appears to be the case regard- 
less of whether the projectile was originally solid and 
remains solid during the collision (low-speed solid-to- 
solid collisions), whether the projectile was originally 
solid but flows during or as a result of the collision 
(high-speed solid-to-solid collisions), or whether the 
projectile was originally a liquid drop (liquid-to-solid 
collisions). At extremely high (meteoric) collision 
velocities it appears that an originally solid target 
may be expected to behave like a liquid dining the 
collision because of the enormous impact pressure 
developed. With regard to collisions of this kind 
Opik [13] has stated that the aerodynamic (1/2 pV 2 ) 
pressure at the penetration of a meteor into rock is 
more than 1,000 times the plastic limit of steel; he 
has hypothesized that all solid materials under such 
pressures must behave like liquids and that the 
problem of meteor impact is the case of the impact 
of a liquid drop of one density into a liquid medium 
of a different density. Under such conditions the 
penetration or pit depth may become independent 
of the collision velocity, or may become essentially 
so, and the pit-depth-versus-velocity curve may ap- 
proach a horizontal line parallel to the velocity axis. 
A trend in this direction may be indicated by the 
flattening of the experimental pit-depth-versus- 
velocity curve for lead. Whenever this condition is 
realized, either partially or completely, the pit- 
depth-versus-velocity equation developed in this 
paper will no longer be applicable. 

3.4. Collision of Copper Spheres Against Copper 
Targets 

The data for impingement of 0.483-cm copper 
spheres against copper targets [11] are shown in 
figure 19 along with the theoretical curve calculated 
by eqs (26) and (27) using physical constants for 
electrolytic copper (table 2). 

The empirical pit-depth-versus-velocity curve for 
copper appears as though it may be flattening at the 
highest collision velocities for which data were ob- 
tained, as was found in the case of lead. If this is 
indeed true, it raises the question as to what deter- 
mines the flattening of the curve. In the case of 
lead it happened at collision velocities equal to the 
speed of sound in lead. If it is happening for copper 
at velocities of 16X10 4 to 18X10 4 cm/sec, it is 
happening at velocities that are about half the speed 
of sound in copper. The speed of sound in infinite 
medium for electrolytic tough pitch copper is 46.91 
X 10 4 cm/sec. It does not, furthermore, seem to be 
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related to the melting point of the target material 
because copper has a high melting point and lead 
has a very low one. A possible explanation is the 
relative susceptibility of the target material to flow 
by translational slip of the target atoms on their 
lattice planes at high rates of loading. With regard 
to the use of copper as a liner for shaped charges, 
Rinehart and Pearson [14] state that copper flows 
readily at high rates of loading. 



:t 


111111 


1111111 




' ' ' 


' 


.6 









o °xf 


- 


.5 










- 


£ 
°-.4 


- 


go, 






- 


Q. 












f\3 


- 








- 


E 












.2 

.1 










_ 







, 


, 







6 8 10 12 

IMPINGEMENT VELOCITY , Cm /sec 



20 X I0 4 



Figure 19. Calculated curve for collisions of 0.483-cm 
copper spheres against copper targets. 

O observed depth for 0.483-cm sphere. 

3.5. Collisions of Zinc Spheres Against Zinc Targets 

No value of the dynamic compressive yield 
strength of zinc was obtainable. One of the possible 
uses of the pit-depth-versus-velocity equation re- 
ported in this paper is the calculation of the dynamic 
compressive yield strength from pit-depth-versus- 
velocity data. By trial it was found that to obtain 
an acceptable value of the intercept velocity for zinc 
with the pit-depth-versus-velocity data of Partridge 
et al. [11], the value of the dynamic compressive 
yield strength for zinc would be 1.546X10 9 d/cm 2 . 
The experimental pit-depth-versus-velocity data for 
zinc and the line calculated from eqs (26) and (27) 
using this value of the dynamic compressive yield 
strength and the physical constants given in table 
2 are shown in figure 20. The agreement between 
the experimental points and the theoretical curve is 
reasonably good. The two points obtained at the 
highest velocities used seem to show a flattening of 
the curve. More experimental data are needed to 
determine whether this is the case. 



Washington, January 29, 1959. 
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Figure 20. Calculated curve for collisions of 0.483-cm zinc 
spheres against zinc targets. 

G observed depth for 0.483-cm sphere. 

3.6. Collisions of Tin Spheres Against Tin Targets 

' Partridge et al. [11] also obtained pit-depth- 
versus-velocity data for collisions of tin spheres 
against tin targets. The tin appeared to crack and 
break off rather than to flow plastically around the 
crater. This interesting fact recalls the similarity 
in behavior of materials at low temperature and at 
high rates of loading. At low temperature tin tends 
to exist in the brittle gray tin form. It is doubtful 
whether the pit-depth-versus-velocity equation can 
be validly applied to tin because of the energy that is 
diverted from pit formation into crack formation 
(see section 2.8). 
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